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G^ I Abstract 
. 

. The phase diagram of non-compact lattice QED in four dimensions with 

I I staggered fermions of charges 1 and —1/2 is investigated. The renormalized 
charges are determined and found to be in agreement with perturbation the- 

• ory. This is an indication that there is no continuum limit with non-vanishing 

^ i renormalized gauge coupling, and that the theory has a validity bound for 



every finite value of the renormalized coupling. The renormahzation group 
flow of the charges is investigated and an estimate for the validity bound as a 
function of the cut-off is obtained. Generalizing this estimate to all fermions in 
the Standard Model, it is found that a cut-off at the Planck scale implies that 
has to be less than 1/80. Due to spontaneous chiral symmetry breaking, 
strongly bound fermion-antifermion composite states are generated. Their 
spectrum is discussed. 



1 Introduction 



Interest in non-perturbative investigations of QED has a long history. From per- 
turbation theory there are indications (Landau pole that the cut-off can only 
be removed from the theory if the renormalized charge vanishes. Also, it is a yet 
unsolved fundamental question [0, why the fine structure constant has the value 
a ~ 1/137. At strong coupling, QED exhibits a chiral symmetry breaking phase 
transition of 2nd order ^, |^, where tightly bound fermion-antifermion pairs are 
generated. There one could expect a deviation from the charge screening behaviour 
known from perturbative QED. The phase transition makes strongly coupled 
gauge theories interesting also for applications in technicolour theories because of 
large anomalous dimensions of the operator . For several years, the possibility 

of a non-trivial continuum limit, or a non-trivial ultraviolet-stable fixed point of the 
Callan-Symanzik /? function, has been investigated in the non-compact formulation 
of QED on the lattice (see for example p|-||l6|). 

Studying non-compact lattice QED with dynamical fermions, some groups find 
non-trivial scaling behaviour |jl2|, |17|, 0, others find their critical exponents to 
indicate triviality [0, 11^ . The QED (5 function has been studied non-perturbatively 



on the lattice |Tq, EO], and using Schwinger-Dyson equations [21 1. It turns out to be 



in agreement with perturbation theory and to show no indications of a non-trivial 
fixed point. If QED is trivial in the limit of infinite cut-off, there is a maximal cut-off 
corresponding to every finite value of the renormalized charge. One is interested in 
the size of this validity bound. A rough estimate of it has been made in a lattice 
study of QED with one charge . 

In nature one finds several differently charged types of fermions. In the presence 
of several charges, new non-perturbative phenomena may arise which could affect 
the phase structure. An important question is whether chiral symmetry breaking 
sets in at the same value of the bare coupling for all fermions or whether some 
fermion species can be massless and others in the chirally broken phase at the same 
time. The phase diagram of such a system has to be investigated with respect to a 
physically interesting continuum limit. In QED with one charge, neutral pointlike 
Goldstone bosons and scalar particles are generated due to the chiral symmetry 
breaking phase transition. One could expect that in the two-charge model there 
is a larger number of neutral pointlike bound states which could in principle carry 
colour charge and have an effect on the /3 function of QCD. Moreover, electrically 
charged bound states may appear which could change the behaviour of the QED /3 
function and push the validity bound towards much lower energies. It is of interest 
to see whether charge renormalization is in agreement with perturbation theory and 
which states give relevant contributions. 

To investigate non-perturbative phenomena in the coupling between different 
species of fermions, a model with two species of staggered fermions with a ratio 
of their charges of —1/2 ('two-charge model') was studied, comparable to u and d 
quarks whose strong and weak interactions are switched off. 
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In section ^ of this paper a description of the model and details of the simulation 
are given. The phase diagram, obtained from the chiral condensates, and the scaling 
behaviour are discussed. The charged sector of the model is presented in section ^. 
The renormalized charges are determined non-perturbatively using current-photon 
correlation functions and are compared with renormalized lattice perturbation the- 
ory. There is good agreement which leads to the conclusion that the model is trivial. 
Using perturbation theory, the renormalization group flow of the charges is deter- 
mined, which leads to an estimate of the validity bound resulting from triviality. 
Spectrum and renormalization group flows of fermion-antifermion composite states 
are discussed in section ^. 

2 Phase diagram and scaling behaviour 
2.1 Action and Simulation Details 

The 'two-charge model' contains a non-compact gauge field with the action 

F^,,{x) = A^A,{x) - A,A^{x), (2) 

and two sets with each four flavours of staggered fermions which couple to the 
gauge fields with couplings 1 and -1/2, corresponding to u and d quarks with only 
electromagnetic interactions: 

= J2 {Xu{x)Mu,xyXu{y) + Xd{x)Md,xyXd{y)} ■ (3) 

x,y 

The coupling (3 is related to the bare electric charges Cfc by = l/{ckek)'^, k = 1,2. 
The lattice spacing a is set to 1. The fermion matrices are given by 

Mk,x,y = Dk,xy + mkSxy 



Dk,xy (•^) 



2 M 



k = u,d ; Cu = I, Cd = -1/2. 

In the limit rrid oo this model goes over into non-compact QED with one 
charge For the gauge fields, periodic boundary conditions in all four di- 

rections were chosen, for the fermions periodic spatial and antiperiodic temporal 
boundary conditions. The simulations were performed on lattices of size 8^ X 12. 
From simulations with one charge, one expects that for the chosen values of j3, 
rrtu and finite size effects are small |r3|. For each simulated point {i3,mu,md) 



0(1000) configurations in equilibrium were generated using a Hybrid Monte Carlo 
algorithm. Every fifth was stored for spectrum and charge calculations. 
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Staggered fermions are a useful choice for studying chiral symmetry properties at 
finite lattice spacing. The action (j^) of the two-charge model has for the kth species 
of fermions a chiral U{l)v x U{1)a symmetry, if = 0. The order parameters are 
the chiral condensates 

(Tfe = ixkXk) = -{TrM^'), (5) 

which are computed using a stochastic estimator [^]. Simulation results for the 
chiral condensates are shown in tables |1| and ^ 



2.2 Determination of the critical points 

In non-compact QED with one set of staggered fermions the chiral condensate a 
is consistent with a mean field like equation of state with logarithmic corrections 
motivated from a linear a model The parameters in the equation of state are 

expanded in a power series in the reduced coupling (1 — /3//3c), where j3c denotes the 
critical coupling. The logarithmic corrections are only expected to become important 
very close to the critical point due to renormalization effects. It is expected that 
they become relevant also here if one goes closer to the critical point. From the 
results in tables |1] and 0, and as illustrated for (3 = 0.18 in figure |I], it appears 
that the chiral condensates are fairly independent of the other fermion's bare mass. 
Thus to determine the critical points in the two-charge model an Ansatz with two 
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Figure 1: Dependence of au on m = rrid (left) and of cr^ on m = (right). Circles 
denote m = 0.02, squares m = 0.04 and diamonds m = 0.09. (3 = 0.18. Statistical 
errors are smaller than the symbols. 

uncoupled mean field like equations of state with logarithmic corrections for each 
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fermion is used: 



111^"= la, 



In lui, ^1 ' 



k = u,d. 



(6) 



For the determination of the renormahzation group flow it is desirable to know the 
chiral condensates as a function of (/?, m„, rrid) in the whole parameter space. It turns 
out to be possible to approximate the chiral condensates for all (3 with equations of 
state (P), using the following expansion of the couplings: 



3 

cuj 1 



(7) 



and 



rf(l-/5//3,,), 



3(1) 



(8) 



(1,3) ^(0,1,3) 



, i/^, and pk are fit parameters. Including all results for simulated 

.05 .10 .15 .20 




Figure 2: Chiral condensates (open symbols: a^, filled symbols: a^) as a function 
of (3 with a mean field fit with logarithmic corrections. Triangles correspond to ruk 
= 0.02, squares to ruk = 0.04, asterisks to ruk = 0.09 and circles to = 0.16, 
k = u,d. The solid lines show the result of the fits and the dashed lines solutions of 
the equations of state with set to 0. 
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at rriu < 0.09, one obtains (3cu = 0.173. Values of the fit parameters and the fit 
errors are hsted in table |^. For a fit of a^, all results with (3 < 0.1 and < 0.09 
have been used, and the critical coupling is f3cd = 0.047. From varying the range 
of the chiral condensates included in the fit, one estimates the error on j3cu and jScd 
to be approximately 0.001, which is larger than the fit error. Without including 
logarithmic corrections, Pcu comes out to be 0.183(1) and Pcd to be 0.049(1) p5| . 
The cubic term in eq. (^ is included to obtain an approximate description of cr^ 
also in the region where j3 ~ f3cd- It has been checked that with a fit in the range 
0.16 < f3 < 0.22, including only linear terms in the reduced coupling, one obtains the 
same result for Pcu within errors. As seen in figure ^ eqs. (|^) give a good description 
for the chiral condensates in the range of couplings (3 > 0.05. Two or three of the 
same symbols lying on top of each other at the (3 values 0.06, 0.16, 0.18, 0.20 and 
0.21 corresponds to simulations at various values of m„ at a fixed value of rrid or 
vice versa. One further notices that in the regions where the d fermion undergoes 
a transition, the u fermion is that far in the broken region that (Xu is practically 
independent of rriu (as well as of m^). 



2.3 Renormalized fermion masses 

The next step in the investigation of the critical behaviour is the determination of 
the renormalized fermion masses or inverse fermionic correlation lengths. Because 
the fermions are charged, their correlation functions are gauge dependent. For the 
calculation of their expectation value, a technique as described in references p6| , |T^ 



IS 

used. First, Landau gauge is fixed by imposing the following gauge fixing condition: 

Y.A,A^ix)=0, (9) 

where denotes the backward derivative on the lattice. An additional gauge-like 
degree of freedom is the invariance of the action under the local transformation 

Xix) e-^'^'^-^^h.ix), 

Ux) e'^^-"(^)xfc(x), (10) 

with ^ 

(y{x) ='^—nf,x^ ,nf, e Z, (11) 

where is the lattice extent in the /i direction. The lattice average of the gauge 
field 

A, = ^Y.M^), v = Y[L„ (12) 

has a nonvanishing expectation value on our (relatively small) ensembles. By shifting 
it by multiples of 47r/L^, such that it is restricted to the interval (— 27r/L^, 27r/L^], 
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Figure 3: u fermion propagators aX (3 = 0.16 and m„ = = 0.02. Crosses 
correspond to positive values of the correlation function, squares to negative values. 
Solid lines denote the fit. 
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Figure 4: d fermion propagators at /3 = 0.16 and = = 0.02. The meaning 
of the symbols is explained above. 
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this additional degree of freedom was fixed. In the two-charge model this interval is 
chosen twice as large as in a model with fermions of charge 1. This is necessary to 
preserve also the boundary conditions of the d fermion, which couples with charge 
— 1/2. is approximately constant over 0(10 — 20) configurations. Following the 
procedure described in |T3[, the set of data samples at each parameter value was 
divided into subsets of 10-20. The correlation functions were averaged over each 
subset and fitted with the free form of a staggered fermion propagator in a constant 
background field B^. For the u fermion, agrees with the expectation value of 
A^, taken over the given subset, for the d fermion with the expectation value of 
—A^j_/2. The fits were performed with the routine MINUIT. The Ansatz gives a 
good description of the data. An example for this is shown in figure |[ Results are 
given in tables ^ and ^. Another indication that the fermion correlation functions 
behave similar to free propagators with a renormalized mass mu^dR, is obtained 
by comparing the simulation data for the chiral condensates with free propagator 
expressions: 

(XkXk) = ttY. ^ ■ I —,k = u,d, (13) 



V Ef, sin + mij^ 



with the lattice momenta 



Pi — IT-i — !)•••) Li, i — 1, . . . , 3, 

TT 

P4 = — (2n4 - 1), n4 = 1, . . . ,L4, 

Figures || and ^ show that for small masses the results agree quite well with eqs. (0) . 
The fermion wave function renormalization constant is 0(1). Since oc (xkXk) 
for small rrikR, eqs. (^ and (|^ imply that near (3 ~ the renormalized u mass 
scales according to 

m„«ln-i/3(m;^)cxmy3, (14) 
and the renormalized d mass according to 

mdRln~P''{m^^) (X nid. (15) 

The scaling behaviour of the d is thus in this region close to the perturbative be- 
haviour, which is very different from the behaviour of the u fermion in this region. 
For very large (3 both renormalized masses follow eq. ([T5|). In the neighbourhood of 
j3 ~ jScd, eqs. (|^) and (|^ indicate that 



rridR 



ln-V3(^-^) oc mT- (16) 



In this region the difference between the bare and renormalized masses of the d 
becomes small. 
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Figure 5: Chiral condensate {XuXu) for /5 > 0.14 plotted against the renormalized 
mass rriuR- The hne denotes eq. (0). 
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Figure 6: Chiral condensate {xdXd) plotted against mdR. The hne denotes eq. (|1^). 
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2.4 Phase Diagram 



Figure |^ shows a sketch of the phase diagram. The agreement of the results with eq. 
(^) suggests that there are two distinct regions of chiral symmetry breaking for both 
fermion species. This does not correspond to the expectations in a confining theory. 
If chiral symmetry breaking in non-compact QED was related to confinement, one 
would expect both species of fermions to develop chiral condensates in the same 
time. The end points of the phase boundary of the u, which separates the symmetric 
phase of the u on the rriu = surface from its broken phase, are given by (3cu and 
f3c, the critical point of QED with one set of u fermions |jl5[ which is the limit of 
the two-charge model if rrid oo. In the presence of more charges, the critical 
point is shifted towards stronger coupling, so jScu is slightly smaller than (3c- In the 
limit rriu —>■ oo, which corresponds to QED with one set of charges with coupling 
P = one expects this to occur at /3c/4 (~ 0.046, using the result of [1^). In 
the two-charge model a value very close to this is obtained, Pcd — 0.047. Below 
Pcu, no continuum limit with two fermion species is possible. For the investigated 
parameter values at /3 < 0.14, the renormalized u masses are 0(1) in lattice units, 
which means the u fermion is in this region practically not present in the spectrum. 
In the region (3 > f3cu both fermion masses can go to zero, so this region is the most 
interesting candidate for a continuum limit of the model. 





Figure 7: Sketch of the phase diagram. 
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3 The renormalized coupling 



3.1 Charge determination on the lattice 

The Ward identities ensure that the charge renormahzation is entirely determined 
by the wave function renormahzation of the photon: 

el,dR = Zsel^a- (17) 

Z3 is given by the zero momentum hmit of the gauge invariant part of the photon 
propagator: 

m = {Mk)A,{-k))i^^^, (18) 

/i,A:|fe^=const, 

Z3 = \iuiDik). (19) 



The sum in eq. (|T8|) runs over aU directions /i and for each /i over aU choices 
of k with fixed fc^ and k^ = 0. Due to the lattice, the momenta enter the photon 
propagator as 

k^ = e'^^ - 1 and k^ = kji^. (20) 

The right hand side of (0) turns out to be strongly fluctuating and inappropriate 
for an extrapolation to /c 0. Thus for the calculation of D{k) a method analogous 
to references |TB|, ^ is used. The photon propagator D{k) is re-expressed using 
the Ward identities p8| in terms of a correlator between the gauge field and the 
fermion current: 



[21] 



k^=0 



where in the two-charge model 



•^'^^^) = E {iXuiy)M^^y,Xu{z) + Xdiy)Md,y,Xdiz)} . (22) 

The correlator in eq. (^Tj) has less fluctuations and could be used for an extrapolation 
of D{k) to A; — > 0. The fermion current was computed using a stochastic estimator 
with 30-75 inversions of the fermion matrices. 



3.2 Comparison with perturbation theory 

Taking contributions of both fermions into account, in one-loop perturbation theory 
the vacuum polarization tensor has the following form: 

1/(3 n^,(A;, m„fi, m^n, V) = el n«(fc, m^^, V) + el n«(A;, m,^, V), (23) 
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where Il^^}{k,mji,V) is the vacuum polarization tensor for a single set of fermions 
with the renormalized mass on a lattice of volume V. Projection onto the gauge 
invariant part yields 



D{k) = l + eln^'\k,m^R,V) + elU^'\k,maR,V), 



(24) 



where U.^^\k,'mii,V) is the one loop vacuum polarization function for one set of 
staggered fermions: 



n(i) {k, mR,V) = j:; [n« {k, ruR, V) - n« (0, rriR, V) 



k^ 



(25) 



The second term on the right hand side occurs because n^^^(0, m^, V") 7^ for a 
finite V. This would correspond to a finite photon mass, and the term is subtracted 
off. Here, a fixed ji has been chosen with k^ — 0. Finally one obtains for the photon 
propagator in the two-charge model: 



P 



^ pR + n(0, rriuR, rridR, 00) - n(A;, uIur, nidR, V). 



m 

Extrapolating this to F — > 00, A; — > 0: 

(3r^ p- n(0, niuR, rUdR, 00), 



(26) 



(27) 



gives the perturbative relation between the bare and the renormalized coupling. 
Combining the last two expressions, one gets the fit formula for the Monte Carlo 
results for D{k): 



(5r + n(0, rUuR, TUdR, 00) - n(A;, tUur, rudR, V). 



(28) 



The renormalized coupling Pr is the only free parameter in the fit. 

For the calculation of n(0, ttIur, uidR, 00), the effect of the zero modes of the 
gauge field has to be taken into account. Their contribution is important where 
fermions are light, as are the d fermions at the simulated parameter values with j3 ~ 
(3cu- To calculate the perturbative vacuum polarization functions, the background 
fields are integrated over: 



n(A;, rUuR, rudR, V) = 



J d'^Adet{Du + rriu) det{Dd + md)J2p p{p, k, A,muR,mdR,V) 
J d'^A det(L'„ + niu) det{Dd + ma) 



(29) 

Also in the presence of constant background fields the fermion determinant of free 
fermions can be written as a product of contributions with definite momentum: 



det Mi = det( A + J^ij) = JJ 

ki 



"^i + sin^(A;j,^ + qA^) 



, i — u,d, 



(30) 
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Figure 8: Photon propagator, (3 = 0.17, niu = rrid = 0.02. The circles describe 
the data and the asterisks the fit with the perturbative formula. The extrapolation 
to — s> is denoted by an asterisk. The variable K2 on the x axis is defined as 
K2 = 576/(27r)2fc2 to be an integer. 



with Cu = I and q = —1/2. Choosing /i = 3 in eq. 
follows: 



where 



and 



the function p looks as 



E 



-2- 



E^^3 (cos(fc,,^) sin(p^/2) 
D{ki - p, miR)D{ki + p, m^^) 

1 



+ sin' (A;,,3-P3)/2 



D{ki-p,min) D{ki+p,miR) 
D{k,mii) = ml + Y,sm\k^/2) 

ku = ku + A and kd = kd - A/2. 



cos^(A;i,3),(31) 



(32) 



(33) 



The integral over the background fields was evaluated using a Monte Carlo method, 
representing the background fields through Gaussian random numbers and calculat- 
ing the weight of the fermion determinant in the presence of these background fields. 
Figure | shows an example for a fit of D{k) with one-loop perturbation theory in 
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Figure 9: Renormalized couplings in the region of small niuR and rridR. The line 
denotes the small mass limit of one-loop perturbation theory. 



the presence of background fields. The oscillations of D{k) are caused by the depen- 
dence of k on the direction due to the asymmetry of the lattice. The renormalized 
couplings are listed in tables ^ and |^. 

The renormalized coupling can be calculated directly in perturbation theory, 
using eq. (^). One would like to compare the data with the perturbative result. 
Since this is of interest especially in regions of a high cut-off, i.e. for small masses, 
the limit muR,mdR -C 1 of the vacuum polarization function (pT]) is taken, which 
gives 

(3r- (3 = -n(0, TUuR, rudR, oo) = - — ln{mlj^mdR) + 5c/4. (34) 

From QED with one species of fermions with charge 1, it is known that c ~ 
0.0210 |T5[. Figure ^ shows that in the small mass limit the renormalized charges 
indeed follow the logarithmic behaviour as in the right hand side of eq. (|3^) . For 
larger masses, 11 is no more a function of m^^mdR alone, but the simulation results 



are still in agreement with perturbation theory |3y]. There does not seem to be an 
effect on charge renormalization from possible charged bound states. The agreement 
with perturbation theory indicates that this model is trivial in the limit of infinite 
cut-off, and that all renormalized charges vanish if only one fermion mass goes to 
zero. The renormalized coupling vanishes in those parts of the phase diagram (see 
figure 0) where or ad equal zero. 

One would now like to find a relation between a finite renormalized coupling 
and a corresponding cut-off. The ratio of the renormalized masses has to be kept 
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constant while the cut-off is varied. Re-expressing 
and exponentiating, one obtains: 



in terms of i? = muR/m^R 



rriuR = 1.365 exp 



6n' 



(35) 



From this equation one will be able to estimate the maximal cut-off belonging to 
each renormalized charge. For this, one needs to know which f3 coordinate the point 
with maximal cut-off on a renormalization group flow line with fixed (3r and R has. 
This flow line will be the intersection between the surface with fixed R and the one 
with fixed Pr in the three dimensional parameter space. 



3.3 Renormalization group flow of the coupling 

It is not feasible to cover all the regions of the phase diagram that are of interest 
with simulations. For a determination of the renormalization group flows of the 
charge one therefore has to make use of extrapolation methods. If one is able to 
find a functional dependence of the renormalized mass on the bare parameters it 
is possible to calculate the renormalized coupling using perturbation theory from 
eq. (^). The chiral condensates can be approximated as functions of jS, rriu and 
nid through the equations of state. Figures ^ and |^ suggest that the renormalized 
masses can be expressed as functions of the chiral condensate, where for small masses 
there is agreement with the free propagator relation eq. ([I3|), and the leading term 
is linear. For the renormalization group flows one is primarily interested in regions 
where both fermion masses are smaller than 1, so m^/j is fltted for (3 > 0.14 with 
the Ansatz 

rriuR = PiiXuXu) + P2{XuXuf ■ (36) 

Here, x^ per degree of freedom is 1.5. The following fit gives a good description for 
the d mass for (3 > 0.14: 

mdR = QiiXdXd), (37) 

with per degree of freedom equal to 2.2. A cubic term is not needed here. The 
fit parameters P and Q are given in table |^. As illustrated in figure one thus 
obtains a good description for the renormalized masses. To obtain the renormalized 
masses as a function of the bare parameters, first the equations of state were solved 
explicitly for the chiral condensates using Miiller's method for calculating zeros 



of non-linear functions. The results were thus inserted into the fit equations (|36D 
and ( p7| ) respectively. 

Using this method, a grid of values for the mass ratio R in the phase diagram is 
generated and, using renormalized perturbation theory, also for f^R. Using a graphics 



package, surfaces of constant Pr and R can be drawn. Figure |Tl] shows a surface Sp^^ 
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Figure 10: Renormalized masses as a function of the chiral condensates. The squares 
denote simulation results and lines the fit according to eqs. ( ^61) and (^), respec- 
tively. 
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(black) with a constant (3^ and its intersection with a surface Sr (grey) of constant 
R of 0(1). 

For large m„ and both renormalized fermion masses are large. Therefore 
[3r ~ /3, which implies that surfaces with fixed [3r are nearly perpendicular to the /3 
axis. The behaviour at small m^, md can be read off from eq. (0). If is lowered, 
the surfaces bend over until they end on the m„ = plane below the line where the 
chiral symmetry of the u fermion is broken. The larger (3r is, the closer the end line 
of the surface on the ttIu = plane comes to the critical line of the u, except for 
very small m^^. In a whole range oi (3 < (3cu and small bare masses, m^^ is much 
smaller than and for constant nearly constant. So surfaces with constant 
(3fi bend at small m^^ towards small f3, as the black surface in figure ^ indicates. 
They intersect the = plane in the broken phase of the d fermion, and thus 
cut the (3 axis at /5 < Pcd- The important point to note is that 5*^^ intersects the 
surfaces rriu = or rrid = only in regions where both renormalized fermion masses 
are non-zero. 




Figure 11: Surfaces with /Sr = 0.27 (black) and R = 5 (grey). The little spheres 
denote Pcu and Pcd- 

Surfaces of constant R end for (3 > f3cu on the (3 axis. The scaling behaviour 
of the masses in eqs. ([T^ ) and ([15|) implies that at /? = /?c« these surfaces obey the 
relation m„ ~ cm^, where c is a constant. If (3 is lowered past (3cui increases 



quickly. As indicated in figure ITl|, the end line of the surfaces with constant R on the 



= plane therefore has to turn away from the P axis towards a finite nid. For an 
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estimate of the validity bound, the intersection hne of both surfaces is of interest, 
as it represents a flow line where both R and (3^ are constant while the cut-off scale 
is changed. It ends on the m„ = plane at /? < [3cu- As one moves down along the 
intersection line, the cut-off becomes larger, but never reaches infinity. 



3.4 Validity bound 

In non-compact QED with four flavours of staggered fermions of charge 1, the (3 
function goes, in the limit of small renormalized mass, over into the prediction of 
one loop continuum perturbation theory for four fermion flavours of charge 1 [|T5| . 



In the two-charge model with four flavours of charge 1 and four flavours of charge 
— 1/2, a corresponding relation is fulfllled (0). This suggests that in the presence 
of species of fermions the renormalized charge can be approximated with the 
formula 

- /5 = - E In (^^iR) + E (38) 

^=l i=l 

In this section the lattice spacing a is kept explicitly, to illustrate the dependence on 
the cut-off. In QED with one charge and in the two-charge model it was reasonable 
to approximate the bare charge corresponding to the point of maximal cut-off by the 
critical coupling of the fermion with the strongest coupling. This critical coupling 
is dependent on the number of dynamical fermions. Comparing jSc, determined 
using mean fleld equations of state or mean fleld equations of state with logarithmic 



corrections, in the quenched case, in models with four flavours of charge 1 [0], with 
eight flavours of charge 1 and the two-charge model, one flnds the following 
behaviour: 

AT 

/3c(iV species) = /9c(quenched) — eE Qi) (39) 

i=l 

where e ~ 0.01. 

In the Standard Model, one has three generations of fermions of charge 1, three 
generations and colours of charge 2/3 and three generations and colours of charge 
1/3. Expressing the electron mass in units of the cut-off and the other fermion 
masses in terms of ratios with the electron mass, one obtains 

am,R > 3.68exp |-^/?«| x (R^Rr)'^' {RuRcRtf {RARbf , (40) 

where R^ = rrieR/m^ji etc. The relation between the lattice spacing a and the cut-off 
A is approximately given by the relation: 

A~l/a. (41) 

Using the physical fermion masses and the physical value of the flne structure con- 
stant 1/137, one gets the cut-off 

A < 10=^2 GeV. (42) 
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This is much larger than the Planck scale (10^^ GeV). However, if there are more 
charged particles, e.g. due to supersymmetry, the exponential dependence on the 
number of charged particles might cause this cut-off to become considerably lower. 

If on the other hand one sets the validity bound of QED to be at the Planck 
scale one obtains the upper bound on the fine structure constant 

an < 1/80, (43) 



which is surprisingly small. One also has to note that in eq. (^31) the effect of 
charged bosons is not yet included. 

4 Composite states 

4.1 Lattice operators and fits 

Correlation functions of scalar and pseudoscalar uu and dd states were investigated: 

CuAt) = (^T.^i(^^t)Xu{x,t)xu{x,t)xu{0)Xu{0)^ , ^ = 1, 2, (44) 
CciAt) = ^*(^' t)Xd{S, t)xd{x, t)xMXd{0)^ , 1 = 1,2. (45) 

The sign factors Sj determine the lattice representation a state belongs to . The 
corresponding continuum quantum numbers and states in QCD terminology are 
listed in table 0. To reduce statistical fluctuations, for each correlation function 32 
local sources distributed over the lattice were used. Masses were determined by fits 
to the formula: 

C{t) = A, (e-^^* + e-^^(^"*)) + A^i-iy (e"^^* + e'^^^^-*)) + A3 + {-IYA,. (46) 

Here, Ei is the energy of the lightest pseudoscalar and E2 the energy of the lightest 
scalar state. The constants A^ and correspond to single fermions which propagate 
in the time direction around the lattice (see |^). Here they were not needed for 
fits of uu states in any parameter region investigated. Masses of neutral composite 
states for /5 ~ Pcu are presented in table || and for /5 ~ Pcd in table ||. 

For correlation functions of type 2, fits were done with A2 set to zero. In this 
channel only the Goldstone pion contributes, since states with the quantum numbers 



cannot be realized in the quark model. Figure shows correlation functions 
of type 2 at /? ~ jScd and figure at /3 ~ Pcu- For uu states a fit interval tmin/tmax 
= 1/11 was chosen. Close to (3cd the mass of the vr^ is nearly independent of the fit 
interval, which indicates that there is a good overlap of the pointlike interpolating 
field with the pion. For f3 < 0.19, the masses of the uu pions are smaller than 2m„/j, 
thus one can speak of bound states, which become lighter the deeper one goes into 
the broken phase. For j3 > 0.20 the vr^ masses lie above 2muR- This is an indication 
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that the pion is not bound any more, instead there is possibly in the infinite volume 
limit a resonance in the pseudoscalar channel. It has to be noted that when a 
spectrum with many states is fitted with a single exponential, the fit result might be 
an average between the lowest lying state around 2mR and the excited states ||33[| . 

For dd states the fit range was at 13 < 0.10 chosen to be tmin/tmax = 1/11. 
For larger good fits could not be obtained unless the constants ^3 and ^4 were 
included. and A4 are about O(10~^) smaller than Ai. The fit interval was in 
this /? range chosen to be tmin/tmax = 2/10. The vr^ masses are around Pcd below 
2mciR, for P > 0.05 there is no bound pseudoscalar dd state. 

Figures |1^ and |15] show typical correlation functions of type 1 for couplings close 
to Pcu and Pcd- For P values much lower than Pcu there is no clear signal for the 
au correlation function. The results shown in tables ^ and ^ were obtained without 
including the pseudoscalar contribution in the fit. Apparently, the pseudoscalar 
state does not give an important contribution to correlation functions of type 1, its 
amplitude is about an order of magnitude smaller than the one of the scalar state 
and per degree of freedom is the same whether the pseudoscalar state is included 
in the fit or not. The difference in the results from both fits is of the same order 
as the error due to the choice of the fit interval. Fit intervals were tmin/tmax = 
2/10. The du energies lie for (3 values between 0.15 and 0.17 slightly below 2muR 
so that they might be bound in this region. For larger f3 they get larger than twice 
the renormalized fermion mass. In the parameter region studied, all ad masses are 
larger than 2mdR- From the simulation results at /3 > 0.075, one observes that far 
in the symmetric phase tt^ and ad states are degenerate, as one would expect from 
restoration of chiral symmetry. 

Correlation functions for charged composite states [ud) have been calculated for 
P values around Pcu, Pcd and for /3 = 0.12. For small P the correlation functions 
suffer from bad noise problems, and the signal appears to fall off with an energy 
larger than the inverse lattice spacing. For large P there is a good signal, but the 
correlation function falls off with an energy which is 0(1) in lattice units and thus 
much larger than the sum of the renormalized fermion masses. At P = 0.12 there is a 
signal in the channel considered and the energy is of the order of rriuR. However, the 
renormalized u mass is large (0(1)) at those values of the coupling and cannot be 
determined reliably because correlators are very noisy. Thus no evidence for charged 
bound states, with small energies in the limit of large cut-off, has been found. 



4.2 Renormalization group flows 

In section ^ it was discussed that in the two-charge model the cut-off cannot be 
removed if the renormalized coupling is kept finite. In such cases there is always the 
question if there are parameter regions where physics can be kept fairly constant 
while the cut-off is changed. In QED with one set of charges, in general flow lines 
of constant mass ratios involving composite states cross flow lines of constant renor- 
malized charge even before the maximal cut-off resulting from triviality is reached. 
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Figure 12: Correlation functions type 2, /9 = 0.05, m„ = m^^ = 0.02, 7r„ (left) and 
71 d (right). Crosses denote positive values of the correlation function, the solid line 
denotes the fit. 



5 10 5 10 




5 10 I 5 10 I 



Figure 13: Correlation functions type 2, /3 = 0.17, m„ = = 0.02, tt^ (left) and 
TTd (right). The symbols have the same meaning as above. 
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Figure 14: Correlation functions type 1, (3 = 0.05, m„ = = 0.02, a^. Crosses 
denote positive and squares negative values of the correlation function, the solid line 
denotes the fit. 
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^ . ^ ^ \ ^ , ^ ^ \ ^ , ,0, I ^ . ^ ^ \ ^ . ^ ^ \ ^ 

c(i)° 
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Figure 15: Correlation functions type 1, f3 = 0.17, m„ = = 0.02, (left) and 
ad (right). The symbols have the same meaning as above. 
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Only in the sector with very small charges and masses the flow lines are nearly paral- 
lel and physics can be kept nearly constant. There are indications that the situation 
can be improved if four-fermion interactions are included In the two-charge 



model the two species of fermions interact with each other only weakly and one 
expects the behaviour of flows belonging to the d to be completely different from 
the behaviour of flow lines belonging to the u. 

The physically most interesting region in the phase diagram is around (3 ~ (3cu 
and at small bare masses. Looking at table it seems reasonable to assume that 
energies of composite states of one fermion only have a very small dependence on the 
bare mass of the other fermion and that thus the flow lines in the plane m = m„ = 
can easily be generalized to the two dimensional flow surface. To obtain flow lines 
in the m„ = md plane, an interpolation between the grid of actual simulation results 
in this plane is performed. For this, the following dependence of mass ratios S on 
the simulation parameters is assumed: 

In 5 = a + h(3, m fixed (47) 
InS = c + (i In m, /3 fixed. (48) 

The Ansatz in eq. ( ^7[ ) is motivated by the logarithmic relation eq. (pSf ) between 
renormalized masses and the charge, whereas eq. (ff3) is motivated by the scaling be- 



haviour expected from the equations of state. Figure |1^ shows lines with a constant 
ratio S = muR/rriT^^ in the rriu = rrid plane. S is varied in steps of 0.1 from 0.4 to 
2.1. The picture suggests that the lines flow into f3cu- Generalizing this to different 
values of m„ and m^, one obtains a picture about the flows as shown in figure |1^. 
Surfaces S = const are expected to end on the m„ = plane at the phase boundary 



of the u. As the thick black line in figure |T7| indicates, one generally cannot keep 
the ratio of fermion and pion masses constant on a line with constant fermion mass 
ratio and renormalized charge, except probably in the perturbative region. 



Figure shows lines with a constant ratio 5* = rriuR/m^^. S is varied in steps 
of 0.025 from 0.3 to 0.6. The lines do not flow into Pcu- One expects that lines of 
constant (3ji and fermion mass ratio will also in general not lie on the surface one 
obtains from generalizing S to the three-dimensional parameter space. However, 
for very small couplings and masses it seems that flow lines with constant coupling 
and fermion mass ratio follow surfaces with constant mass ratios closely and renor- 
malizability is essentially restored. Lines with constant S = mfiR/nij^^ are shown in 
figure |l^. 5* is varied in steps of 0.0125 from 0.025 to 0.1. The lines show that this 
mass ratio is for large couplings fairly independent of irid- 
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,15 .16 .17 .18 .19 .20 .21 

Figure 16: Lines with constant m„fl/m^^ in the = ma plane. The black dot 
denotes (3cu- 




Figure 17: Sketch of a surface with constant and a surface with constant 
muR/rriT^^. The thick black line denotes a constant fermion mass ratio on the surface 
with constant 
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Figure 18: Lines with constant muR/rria^ in the 
denotes Pcu- 



— rrid plane. The black dot 




Figure 19: Lines with constant mdR/rriT,^ in the m„ = rud plane. The black dot 
denotes (5cu- 
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5 Conclusions 



In this paper a lattice study of non-compact QED with two sets of staggered fermions 
with charges 1 (u) and —1/2 (d) ('two-charge model'), is presented. The phase dia- 
gram is obtained from the chiral condensates. They can be described by a fit with 
equations of state of an 0(2) symmetric linear sigma model with logarithmic cor- 
rections to the mean-field equations. Chiral symmetry breaking occurs at different 
values of the bare coupling for both fermions, for the u fermion at /3c« = 0.173(1) 
and for the d fermion at (3cd = 0.047(1). The most interesting candidate for a con- 
tinuum limit of the model is at j3cu, with m,u = rrid = 0. This is the end point of 
the line on the /3 axis where renormalized masses of both fermions are zero in units 
of the cut-off. There are indications that for smaller P the renormahzed d mass can 
go to zero, while the renormalized u mass is finite. If (3 is lowered past Pcd, both 
renormalizd masses are always finite. 

The renormalized coupling has been determined and found to be compatible 
with perturbation theory. Other effects to the charge renormalization like possible 
charged bound states seem not to give a noticeable contribution. The agreement with 
perturbation theory indicates that the renormalized charge of all fermions vanishes 
even if only one becomes massless. An estimate for the validity bound of the two- 
charge model was obtained and generalized to all charged fermions in the Standard 
Model. Including all known charged fermions one gets an upper bound of < 1/80 
if one assumes QED to be vahd up to the Planck scale. 

Study of composite states {uu and dd) has shown that in the neighbourhood 
of the physically interesting point Pcu only the uu states are bound. Masses of dd 
states are of 0(1) in lattice units in this region of the phase diagram. It appears 
that due to the shape of renormalization group fiows of mass ratios one cannot keep 
physics constant even approximately in the investigated parameter region. The 
theory seems to become inconsistent already at scales which are lower than the 
cut-off due to triviality. However there are indications that renormalizability is 
approximately restored in the perturbative region. The situation that the theory is 
in general not renormalizable may be improved by including other operators into 
the action. 
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Tables 



3 


lllu 


"'(/ 


< XuXa > 


< XdXd > 


0.00 


0.09 


0.09 


0.6312(11) 


0.6323(11) 


0.04 


0.04 


0.04 


0.6346(17) 


0.4036(13) 


0.05 
0.05 
0.05 


0.02 
0.04 
0.09 


0.02 
0.04 
0.09 


0.6302(23) 
0.6295(17) 
0.6258(11) 


0.1881(15) 
0.2515(10) 
0.3339(8) 


0.06 

0.06 


0.02 

0.04 


0.04 
0.04 


0.6180(24) 

0.6243(17) 


0.1392(6) 

0.1397(5) 


0.075 
0.075 
0.075 


0.02 
0.04 
0.09 


0.02 
0.04 
0.09 


0.6067(23) 
0.6073(16) 
0.6061(11) 


0.03988(10) 
0.07898(22) 
0.1662(3) 


0.10 

0.10 
0.10 


0.02 

0.04 
0.09 


0.02 
0.04 
0.09 


0.5569(22) 

0.5640(16) 
0.5664(10) 


0.02602(4) 

0.05380(9) 
0.1148(1) 


0.12 


0.02 


0.04 


0.4970(19) 


0.04357(5) 



Table 1: Chiral condensates on the 8^ x 12 lattice, region /? ~ Pcd- 
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/3 


rriu 




^ XuXu ^ 


< XdXd > 


0.14 


0.02 


0.02 


0.4222(19) 


0.01951(3) 


0.15 


0.02 


0.02 


0.3741(18) 


0.01880(2) 


0.15 


0.04 


0.04 


0.4031(12) 


0.03810(5) 


0.15 


0.09 


0.09 


0.4469(09) 


0.08451(9) 


0.16 


0.02 


0.02 


0.3236(17) 


0.01813(2) 


0.16 


0.02 


0.04 


0.3316(16) 


0.03595(3) 


0.16 


0.04 


0.02 


0.3656(12) 


0.01815(2) 


0.16 


0.04 


0.04 


0.3629(12) 


0.03611(3) 


0.16 


0.09 


0.09 


0.4168(8) 


0.08146(8) 


0.17 


0.02 


0.02 


0.2742(15) 


0.01772(2) 


0.17 


0.04 


0.04 


0.3238(12) 


0.03539(4) 


0.17 


0.09 


0.09 


0.3908(8) 


0.07871(7) 


0.18 


0.02 


0.02 


0.2167(14) 


0.01741(2) 


0.18 


0.02 


0.04 


0.2227(15) 


0.03435(3) 


0.18 


0.02 


0.09 


0.2255(14) 


0.07602(7) 


0.18 


0.04 


0.02 


0.2850(11) 


0.01725(2) 


0.18 


0.04 


0.04 


0.2854(11) 


0.03524(5) 


0.18 


0.04 


0.09 


0.2844(12) 


0.07702(8) 


0.18 


0.09 


0.02 


0.3611(8) 


0.01737(2) 


0.18 


0.09 


0.04 


0.3618(8) 


0.03473(3) 


0.18 


0.09 


0.09 


0.3621(8) 


0.07716(7) 


0.19 


0.02 


0.02 


0.1870(15) 


0.01669(1) 


0.19 


0.04 


0.04 


0.2444(10) 


0.03370(4) 


0.19 


0.09 


0.09 


0.3372(7) 


0.07502(6) 


0.20 


0.02 


0.02 


0.1416(10) 


0.01649(1) 


0.20 


0.02 


0.16 


0.1462(15) 


0.1263(1) 


0.20 


0.04 


0.04 


0.2167(9) 


0.03283(3) 


0.20 


0.09 


0.09 


0.3108(8) 


0.07420(6) 


0.20 


0.16 


0.02 


0.3785(6) 


0.01677(1) 


0.21 


0.02 


0.02 


0.1143(7) 


0.01624(1) 


0.21 


0.02 


0.16 


0.1072(7) 


0.1254(1) 


0.21 


0.04 


0.04 


0.1838(8) 


0.03331(6) 


0.21 


0.09 


0.09 


0.2882(7) 


0.07280(6) 


0.21 


0.16 


0.02 


0.3605(5) 


0.01660(1) 


0.22 


0.04 


0.04 


0.1601(7) 


0.03222(3) 


0.22 


0.09 


0.09 


0.2693(6) 


0.07162(6) 



Table 2: Chiral condensates on the 8^^ x 12 lattice, region /3 ~ /3< 
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Pen 


0.1729(2) 


I3cd 


0.04675(7) 


^(1) 


-1.78(1) 




-1.469(9) 


'u 


0.48(3) 


Of 




0(0) 


2.61(2) 


2.41(2) 


0(1) 

u 


-1.3(1) 


0(1) 


-1.42(2) 


0(3) 


-5.7(6) 






Pu 


0.588(6) 


Pd 


0.074(3) 



Table 3: Fit parameters for the equation of state with logarithmic corrections. 



p 


m„ 


rrid 




rridR 


Pr 


0.05 
0.05 

0.05 


0.02 
0.04 

0.09 


0.02 
0.04 

0.09 


5.000(675) 
4.935(143) 
5.616(601) 


0.326(7) 
0.446(4) 
0.611(7) 


0.073(1) 
0.069(1) 

0.063(1) 


0.06 
0.06 


0.02 
0.04 


0.04 
0.04 


2.872(393) 
2.573(382) 


0.258(4) 
0.252(3) 


0.087(1) 
0.088(1) 


0.075 
0.075 
0.075 


0.02 
0.04 
0.09 


0.02 
0.04 
0.09 
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0.085(1) 
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0.10 
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0.04 
0.09 
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0.09 
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0.050(1) 
0.090(1) 
0.209(2) 


0.163(2) 
0.154(1) 
0.139(1) 


0.12 


0.02 


0.04 


1.124(37) 


0.072(1) 





Table 4: Renormalized fermion masses and coupling, region /? ~ Pcd- 
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0.17 


0.04 


0.04 


0.607 


;§) 


0.061( 


1) 


0.274 


;3) 


0.17 


0.09 


0.09 


0.815 


;5) 


0.139( 


1) 


0.250 


[2) 


0.18 


0.02 


0.02 


0.366 
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1) 
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1) 
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'3) 
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1) 
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1) 
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1) 
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[2) 


0.19 
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1) 
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[3) 
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[2) 


0.20 


0.02 
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1) 
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;3) 
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1) 
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1) 


0.342 


;3) 
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0.09 


0.596 
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1) 


0.299 


[2) 


0.20 


0.16 


0.02 


0.807 
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1) 


0.303 


[2) 


0.21 


0.02 


0.02 
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;4) 
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1) 
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[2) 


0.21 


0.02 


0.16 


0.171 


[2) 


0.213( 


1) 
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[3) 
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0.311 


[3) 


0.055( 


1) 


0.367 


[3) 


0.21 


0.09 


0.09 


0.528 




0.125( 


1) 


0.315 


[2) 


0.21 


0.16 


0.02 


0.735 


y) 


0.029( 


1) 


0.318 


[2) 


0.22 


0.04 


0.04 


0.268 


;4) 


0.054( 


1) 


0.386 


;3) 


0.22 


0.09 


0.09 


0.491 


;3) 


0.124( 


1) 


0.328 


[2) 



Table 5: Renormalized fermion masses and coupling, region (3 ^ (3^ 
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Pi 


1.54(1) 


Qi 


1.726(3) 


P2 


3.5(1) 


Q2 





Table 6: Fit parameters for the renormalized fermion masses. 
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continuum states 
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or 


TT 



Table 7: Sign factors for the meson operators. The a particle is provided with a hat 
to avoid confusion with the chiral condensates introduced earlier. 















0.05 


0.02 


0.02 


0.305(1) 


0.431(1) 


0.714(12) 


0.05 


0.04 


0.04 


0.431(0) 


0.556(1) 


0.950(21) 


0.05 


0.09 


0.09 


0.648(0) 


0.783(1) 


1.219(16) 
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0.02 


0.04 
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0.695(2) 
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0.311(1) 
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0.04 


0.04 


0.438(1) 
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0.075 
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0.09 
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0.10 


0.02 


0.02 


0.319(1) 


1.124(3) 


1.126(3) 


0.10 


0.04 


0.04 


0.450(1) 


1.042(3) 


1.052(3) 


0.10 


0.09 


0.09 


0.671(1) 


1.142(3) 


1.173(3) 



Table 8: Energies of neutral composite states, region (3 ~ (3cd- 
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Table 9: Energies of neutral composite states, region /3 ~ /S, 
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